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Quench field sensitivity of two-particle correlation in a Hubbard model 
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Short-range interaction can give rise to particle pairing with a short-range correlation, which may 
be destroyed in the presence of an external field. We study the transition between correlated and 
uncorrelated particle states in the framework of one-dimensional Hubbard model driven by a field. 
We show that the long time-scale transfer rate from an initial correlated state to final uncorrelated 
particle states is sensitive to the quench field strength and exhibits a periodic behavior. This 
process involves an irreversible energy transfer from the field to particles, leading to a quantum 
electrothermal effect. 
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I. INTRODUCTION 

The existence of short-range interaction can induce 
many exotic physical phenomena, and therefore boost 
a theoretical interest in a variety of correlated quantum 
systems. Owing to the rapid advance of experimental 
techniques, ultracold quantum gases trapped in optical 
lattices have provided in the past decade a route to sim¬ 
ulate the physics of different kinds of correlated quan¬ 
tum systems [lj. For few interacting particles, exper¬ 
iments with ultracold atoms have so far demonstrated 
the existence of bound pair (BP) states Q and corre¬ 
lated tunneling phenomena [1]. This stimulate many ex¬ 
perimental and theoretical investigations aiming at the 
correlated particles, which lead to a series of intriguing 
novel phenomena involving BP formation M , detection 
Et dynamics HEMzl , and BP condensate (l8| . 

On the other hand, recent experiments on optical 
lattices [H| have renewed the interest for investigat¬ 
ing the behavior of correlated quantum system after a 
sudden change of Hamiltonian parameters, the so-called 
quantum quench. The corresponding spectacular ex¬ 
perimental results have triggered an intensive research 
on quantum quenches in various systems involving one- 
dimensional Bose systems [20j, Luttinger liquids [2l[ and 
others [22} . Among many remarkable quench dynami¬ 
cal phenomena, it has been possible to observe collapse 
and revival of matter waves with Bose-Einstein conden¬ 
sates in optical lattices [23} . coherent quench dynamics 
of fermionic atoms [24} , as well as non-thermal behaviors 
in near-integrable experimental regimes [25} . A straight¬ 
forward method to realize such quenched systems is the 
instantaneous change of a global or local parameter of the 
system like an external field or the interaction strength. 
In particular, when an external field is applied to the cor¬ 
related few-body system, preservation of the correlation 
between the two particles in matter-wave transport can 
exhibit unusual phenomena that are generally inaccessi¬ 
ble in solid-state systems, such as frequency doubling of 
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Bloch oscillation (BO) of two correlated particles 
dynamic localization [30(, coherent destruction of tunnel¬ 
ing 31] and the fractional BO of atomic pairs [32} . How¬ 
ever, if the external field is applied in a suitable manner, 
tunnelling between Bloch bands becomes possible, which 
leads to a sudden death of the BO (Hi or Bloch-Zener 
oscillation 


In previous work [33} . the dynamical behavior of BP 
has been investigated for the extended Hubbard model 
driven by an external field. It has been shown that the 
nearest-neighbor interaction can spoil the completeness 
of the bound band, which can induce a sudden death of 
the oscillation of BP and break the correlation between 
the particles when the external field is switched on. How¬ 
ever, the underlying mechanism of such dynamical be¬ 
haviors is not clear. Motivated by this question, here, 
we theoretically investigate the underlying mechanism of 
the dynamics of the two-particle correlation in the frame 
work of the Hubbard model driven by a field, and show 
that, the presence of the external field can bring about 
an avoided crossing between the correlated energy level 
and uncorrelated scattering band leading to the correla¬ 
tion and energy transfer among involved states. In this 
case, correlated particles can be pumped into the scat¬ 
tering regime accompanied with the energy transfer from 
the field to particles, rather than preserving the correla¬ 
tion of the particles in the unmixed region in which the 
correlated and uncorrelated state are not mixed. Intu¬ 
itively, when an external field is applied, there can ex¬ 
ist a threshold below which the field cannot induce the 
destruction of the short-range correlation between the 
particles. Furthermore, as the magnitude of the field 
strength is beyond the threshold, the transfer rate from 
correlated state to uncorrelated states gets larger with 
increase of the quench field strength. However, in this pa¬ 
per, we show that the long time-scale transfer rate from 
an initial correlated state to final uncorrelated particle 
states exhibits a periodic behavior with the increase of 
the quench field strength, which can occur in a small 
range of the field. The essential physics of such coun¬ 
terintuitive dynamics is the periodic occurrence of the 
field-induced multi avoided crossings between the corre¬ 
lated and uncorrelated energy levels. 
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FIG. 1. (Color online) Two-particle energy spectrum of the 
N = 3 system as a function of parameter F. The other pa¬ 
rameter values of the system are U = V = —6, and k = 0.4. 
The avoided crossing involves the top two energy levels de¬ 
noted by the red and blue line, respectively, which occurs in 
the vicinity of F = U/2 = —3. 


This paper is organized as follows. In Sec. HHwe intro¬ 
duce the model and discuss its dynamics in the avoided 
crossing situation through a simple case. In Sec. EU an 
analysis of the two-particle spectrum is presented, high¬ 
lighting the existence of multi avoided crossings. Sec. Hvl 
is devoted to investigate the quench dynamics of two- 
particle correlation. Finally, in Sec. |V]the main conclu¬ 
sions are outlined. 


II. MODEL HAMILTONIAN AND THE 
MECHANISM 

For the sake of generality and simplicity, we start from 
a driven one dimensional extended Hubbard model. It is 
worth pointing out that the results shown in the follow¬ 
ing are also applicable to the other form of short-range 
interaction. The model concerned can be employed to 
describe the dynamics of interacting particles on a one- 
dimensional tight-binding lattice driven by an external 
field, which can be defined by the Hamiltonian 

H = Ho + F^jnj, ( 1 ) 

3 = 1 

where the second term describes the linear external 
field served as a quench parameter. And Hq is the 
one-dimensional Hamiltonian for the extended Hubbard 


FIG. 2. (Color online) Plot of the relative distance r as a 
function of parameter F. The red and blue lines denote the 
top two energy levels of the system shown in Fig. [I] It can 
be observed that the two corresponding states exchange the 
correlation through an avoided crossing region. 

model on an N- site lattice 

H q = -kY^ (a]a j+ 1 + H.c.) + ^ n i ~ ^ 

3 =1 3 =1 

+ V J2 n 3 n 3+ 1 > ( 2 ) 

3 = 1 

where aj (a^) is the creation (annihilation) operator for 

boson at site j and rij = is the number operator. 
F is magnitude of an external field. The tunneling am¬ 
plitude and on-site interaction strength are denoted by 
n and U, respectively. And V accounts for the nearest- 
neighbor (NN) interaction. The Hamiltonian © can also 
describe ultracold atoms or molecules with magnetic or 
electric dipole-dipole interactions in optical lattices [40j |. 
In the absence of the external field, the short-range in¬ 
teraction denoted by U and V can induce the particle 
pairing with short-range correlation no matter what mag¬ 
nitude of the particle-particle interaction is [33]. We 
consider a correlated particles state. The spectrum of 
the system after a held quench is dramatically changed, 
which may destroy the correlation between the particles. 
We will demonstrate this point through a simple case. 

We first focus on 3-site chain, which possess a simple 
energy structure but can reveal the underlying mecha¬ 
nism of the field-induced destruction of the correlation 
between the particles. The energy levels of the con¬ 
cerned system as a function of F with the parameter 
U = V = —6 is plotted in Fig. [T] It is shown that there 
are six energy levels. And we focus on the top two energy 
levels, which undergo an avoided crossing at the point of 
U = V = 2 F. When the parameters satisfy the condi¬ 
tion | F — U/2\ <C U = V and F, U k, the dynamics 
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of the correlated particles involving top two energy levels 
is governed by the following effective Hamiltonian 


-Heff = 


\[2 k 2 
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( 3 ) 


which is obtained by the canonical transformation fill ]. 
And |p) = a{ 2 /-\/2 |vac) and |up) = ojaj |vac) denote pair 
and impair state, respectively. In order to investigate the 
two-particle correlation of given state, we introduce the 
average distance between the two particles defined as 


r = y^,r(^l rijm+r |T), (4) 

i,r 


which can measure the correlation between the two parti¬ 
cles. We plot the average distance r as a function of F in 
Fig. [2] in which the blue and red line denote the correla¬ 
tion of top two energy levels, respectively. It can be seen 
from that when the external field F varies to the vicin¬ 
ity of the point U = V = 2F, the two states are mixed 
associated with the exchange of correlation through the 
avoided crossing. Based on this, one can assume that 
the occurrence of the avoided region can lead the corre¬ 
lated particles to displaying a distinct dynamical behav¬ 
ior for the system with and without avoided crossing. To 
demonstrate this point, we consider the time evolution 
of the initial state |T (0)} = |p), which can be readily 
obtained as 


IT (t)) = e lCot {[cos (cot) + i cos 6 sin (cot)] |p) 
— [i sin 0 sin (cot)] |up}} , 


where 
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( 5 ) 

( 6 ) 

( 7 ) 

( 8 ) 
( 9 ) 


Here we focus on how much probability of the evolved 
state |T (t)) will remain in the state |p) as time goes on. 
To this aim, we define the transfer rate from the initial 
paired state to the final unpaired state as 


P (t) = 1 — |(p |T (f)}| 2 = sin 2 #sin 2 (oot ). (10) 


For the case of F = U/2 that the parameter is in 
the avoided crossing region, we have Ci = — 4« 2 /3f7, 
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FIG. 3. (Color online) Evolution of the transfer probability 
P(t). The other system parameters are the same as that in 
Fig. m The initial state is prepared in the state |p). The 
blue and red lines represent evolutions of P ( t ) with the two 
typical values of F = —3, —1, which is corresponding to the 
system with and without the avoided crossing, respectively. 


C 0 = AF - 8k 2 /3U, w = V76k 2 /3\U\, and tan# = 3y/2, 
which lead to P (t) = 18 sin 2 (y/l&K 2 t/3 |t/|) /19. The 
transfer rate as a function of time is plotted in Fig. [3] 
It can be shown that the transition between paired and 
unpaired state occurs with period t r = 3 \U\ 7t/v / 76k 2 , 
which is inverse proportion to the energy gap between 
the involved two energy levels. In their essence, the mix¬ 
ing of the pair and unpair state leads to the exchange of 
their correlation in the region of avoided crossing. On the 
other hand, when the parameter is away from the avoided 
crossing region satisfying the condition of |F — U /2| k, 
we have Ci = U — 2F = 2ui, Co = U/2 + 3 F, sin0 = 
tan# = Cv / 2k 2 /Ci (F 2 — FU), which yield P (t) = 

\Us/2k 2 /C x (F 2 - FU)] 2 sin 2 [(U - 2F) t/2] » 0 owing 
to the condition of |F — U/2] k. This indicates that 
the evolved state will remain in the paired state rather 
than tunneling to the unpaired state as shown in Fig. [3] 
Through the above analysis, it is clear that the existence 
of the avoided crossing is crucial for the coherent motion 
of the correlated particles. However, the structure of the 
avoided crossing for large N system will be something 
different comparing with that of the above N = 3 case. 
In the following section, we will investigate the large N 
system. 


III. ENERGY SPECTRUM AND MULTI 
AVOIDED CROSSINGS 

To investigate the influence of an external field on the 
energy band structure, we first give a two-particle solu¬ 
tion of the Hamiltonian FI q. For the sake of simplicity, 
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we restrict the analysis to the system with U = V < 0, 
which ensures the Bloch band for BP locates below the 
continuum of scattering states. As in Refs. US HI, a 
state in the two-particle Hilbert space, can be expanded 
in the basis set {, r = 0,1,2,...}, with 

K> = E eiK ° ( a 3) i vac > ’ ( u ) 

\tf) = ie iKr/ 2 E e iKja Wi+r l vac ) - ( 12 ) 

v j 

where |vac) is the vacuum state for the boson operator ai. 
Here K denotes the center momentum of two particles, 
and r is the distance between the two particles. Due 
to the translational symmetry of the present system, we 
have the following equivalent Hamiltonian: 

H* = -JkW 2 <0f | + E W) <</f+ 1 | + H.c.) 

i=i 

+u( |^X^| + kf)^f|) ( 13 ) 

in each invariant subspace indexed by K. In its present 
form, H ^ are formally analogous to the tight-binding 
model describing a single-particle dynamics in a semi¬ 
infinite chain with the A'-dependnet hopping integral 
Jk — 2k cos (K/2) in thermodynamic limit N —> oo. In 
this paper, we are interested in the BP states, which cor¬ 
responds to the bound state solution of the single-particle 
Schrodinger equation 


H e q I= fiC I i>K) ■ (14) 

For a given Jk, the Hamiltonian H ^ possesses one or two 
types of bound states Q, which are denoted as \4>k) an£ l 
respectively. Here the Bethe-ansatz wave functions 
have the form 

kl) = C* |^> + (±l) r Cfe-^ \tf) , (15) 

with P > 0. For two such bound states |the 
Schrodinger equation in Eq. m admits 

± e 3/3 + 2'iixe 2 ^ ± ( u 2 k — l) + uk = 0, (16) 


FIG. 4. (Color online) Energy spectrum of the two correlated 
particles as a function of the external held F for the system 
with U = V = —6.24, k = 1. The blue and red region in¬ 
dicates correlated energy level and multi uncorrelated energy 
levels, respectively. It can be observed that a slight variation 
of the external held can give rise to the periodic occurrence 
of the avoided crossing region, in which the correlated state 
and multi uncorrelated states are mixed together. The black 
dashed rectangle denotes a typical avoided crossing, (b) is 
detail with enlarged scale of black dashed marker in (a), (c) 
Schematic illustration of the ideal avoided crossing regions. 
Each of these contains three scattering and two correlated 
energy levels. It should be noted that the slopes of the two 
adjacent correlated lines are slightly different, which is esti¬ 
mated about F. This conclusion can also be applied to the 
adjacent uncorrelated bands composed by multi scattering en¬ 
ergy levels. 


where uk = UI Jk is the reduced interaction strength. 
The corresponding bound-state energy of 1 1 /^) can be 
expressed as 

= ± Jk cosh/3. (17) 

As is shown in Ref. |.33), if \U/n\ is larger than 6, the 
system comprises two complete bands, one of which cor¬ 
responds to two particles bound states. Otherwise, it 
possesses the incomplete bound band associated with the 
vanishing of the gap between the scattering band and 
bound band. Here we want to stress that the preserva¬ 
tion of the correlation between the particles in transport 
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does not depend on the completeness of bound band of 
Hq as the external field is switched on. The dynamics 
of the two correlated particles is determined by the ex¬ 
istence of the avoided crossings of H , which can destroy 
the pair correlation. We will demonstrate this point in 
the following section. 

Now, we switch gears to the spectrum of the Hamilto¬ 
nian H. For the case of U F, k, the energy spectrum 
of the two correlated particles becomes discrete and turns 
out to be composed by two separated Wannier-Stark lad¬ 
ders (39|. In this case, the BP undergoes a perfect BO 
with period is = 2 n/F and preserves the pair correla¬ 
tion, which is analogous to the case of a single particle 
motion subjected to an external field [29, [39j. On the 
other hand, the presence of the external held can not 
only widen the spectrum of the scattering and bound 
band of Hq but also submerge the bound band into the 
scattering band. In this case, surprisingly, we find that 
the correlated energy level can overlap with the uncor¬ 
related band formed by multi uncorrelated energy levels. 
This can be seen from Fig. IHa), in which the blue corre¬ 
lated energy level overlaps with the first red uncorrelated 
band, then in turn encounters the different uncorrelated 
band forming multi avoided regions. In these regions, 
the corresponding wave functions of two types of energy 
levels are mixed accompanied with the correlation and 
energy exchange among the involved states, which can 
prevent the correlated transport of the two particles. It 
is worth pointing out that the existence of multi avoided 
crossing regions is sensitive to the external field F and 
presents a periodic behavior with period F& = 0.0015 in 
such a small range of the external field as shown in Fig. 
Oja) and (c). For the sake of clarity and simplicity, we 
draw schematically the ideal avoided crossing regions be¬ 
tween correlated energy level and scattering band and 
denote the period of the occurrence of these regions in 
Fig. He). One should be noted that the slope of the two 
adjacent blue correlated lines are different and so are the 
red uncorrelated bands, which cannot affect the periodic 
occurrence of the avoided crossing regions in such a small 
range of the external field F. It is presumable that the 
dramatic change of the energy spectrum from avoided re¬ 
gion to unmixed region induced by a slight variation of F 
can be characterized by the dynamics of correlated parti¬ 
cles and exhibits some peculiar behaviors, which will be 
investigated in the following section. 


IV. QUENCH DYNAMICS OF TWO-PARTICLE 
CORRELATION 


To gain insights into the transition of the spectrum 
induced by the external field. We perform a numerical 
simulation to investigate the dynamical behaviors of cor¬ 
related particles. The initial correlated particles state we 



FIG. 5. (Color online) Schematics of the initial correlated 
wave packet projected in the energy space of Ho- The blue 
Gaussian wave packet is centered around Ko = —0.97T of the 
upper bound band denoted by the red line. 


concerned is placed on the upper bound band of FI 0 , 


|T(0)} = A^exp 
K 


(K-Kpf 
2 a 2 


iN A K 


I i>K> , 


(is) 

where A is the normalization factor, Kq and N A denote 
the central momentum and position of the initial wave 
packet, respectively. Throughout this paper, we set the 
parameter K 0 = — 0.97T, a = 0.2, and N A = 36. For 
the sake of clarity, we schematically illustrate the initial 
state IT (0)) in Fig. [5] We consider a quantum quench 
in which the external field strength is suddenly changed 
from F = 0 for f = 0 to F / 0 for t > 0. After the 
quench, the initial state |T (0)) evolves according to the 
new Hamiltonian |T ( t )) = e~ lHt |T (0)). In general, in 
presence of the external field, one would think that there 
can exist a threshold beyond which the initial correlated 
particles state |T (0)) should be pumped into the scatter¬ 
ing band and lose their pair correlation. The correspond¬ 
ing transfer rate from the initial correlated particles state 
to final uncorrelated states gets larger with increase of the 
field strength. However, after a field quench, the dynam¬ 
ics of the correlated particles state will exhibit a periodic 
behavior characterized by two distinct behaviors, owing 
to the dramatic change of the spectrum from the multi 
avoided crossings to unmixed region in a small range of 
the external field. This point will be presented in the 
following. 

Before starting the investigation of the dynamics of 
the two correlated particles state, we would like to study 
the distribution of |T (0)) in the energy space for the 
quenched systems with and without multi avoided cross¬ 
ing regions. The corresponding spectrums of the con¬ 
cerned two systems are plotted in Fig. H For the 
quenched system with multi avoided crossings, the distri- 
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FIG. 6. (Color online) The main distribution of the initial 
BP wave packet in the energy space for the systems with pa¬ 
rameters (a) F = —0.097120, and (b) F = —0.097815, respec¬ 
tively. The parameters U = V = —6.24, n = 1 and N = 111 
are the same for both systems. It is should be noted that we 
restrict the total probability of the initial state distributed 
on the energy space to be 0.95, which can characterize the 
main distribution of the initial state. In (a) the profile of 
energy distribution of the initial state exhibits Gaussian-like 
form and involves less energy, (b) denotes the distribution 
in the system with multi avoided crossing regions involving 
more energy levels comparing with (a). The interval between 
the two adjacent energies related to either correlated states 
or multi scattering states is about F, which agrees with our 
prediction. 


bution of initial correlated state (0)) mainly concen¬ 
trates on the mixed region and therefore refers to many 
energy levels. This is due to the fact that the correla¬ 
tion exchange takes place among the involved correlated 
state and uncorrelated band through the avoided cross¬ 
ing. On the contrary, when the sudden change of the 
field strength F is far away from the avoided crossing 


region, the wave functions of two types of states are not 
mixed, which suppress the tunneling from the correlated 
to uncorrelated states. This leads to the fact that the cor¬ 
responding eigenstates remain their correlation. There¬ 
fore, in contrast to the case of the quenched system with 
multi avoided crossing regions, the initial state 1 1 T r (0)) 
is mainly distributed on the correlated states and con¬ 
tains less energy levels, which are equally spaced with 
an approximately constant spacing given by F. As a 
comparison, we show in Fig. [G] the main distribution of 
the initial correlated wave packet | T (0)) in the energy 
space of H for two typical cases. It can be observed that 
the energy distribution related to the quenched system 
with multi avoided crossing regions involves more energy 
levels comparing with that of the system with the un¬ 
mixed region, which is in agreement with our analysis. 
This transition of the distribution can not only reflect the 
change of the energy spectrum but also result in the en¬ 
tirely different dynamical behaviors of two such quenched 
systems. Thus it can be presumable that the dynam¬ 
ics of the evolved state | T (t)) refered to the quenched 
system without multi avoided crossings will exhibit BO 
with period ts = 2n/F and hold its correlation, which 
is analogous to a single particle motion subjected to an 
external field. On the other hand, for the evolved state 
| U/ (f)) govern by the quenched system with multi avoided 
crossings, the initial correlated particles will be disasso¬ 
ciated and lose the correlation between the two particles 
through the avoided crossing. 

In order to characterize both the above two typical 
behaviors, as in Sec. El we introduce transfer rate T ( t ) 
to measure how much probability of the evolved state will 
remain in the bound states 

T{t) = Y,\(r K \e- iHt \*m \ 2 ■ as) 

cr,K 

We plot Fig. [7]to verify and demonstrate the above anal¬ 
ysis. For the evolved state mainly placed on the multi 
avoided crossing regions, we find that T ( t ) will tend to 
be 0 accompanied by the irreversible increase of r ( t ) as 
time goes on. This indicates that the initial two corre¬ 
lated particles lose their correlation, which can be seen 
from the Fig. Eb). On the contrary, in Fig. [7])a), the 
correlation of the two particles remains strong and the 
corresponding T (t) tends to be a steady value 0.93 in a 
long time scale, which accords with our analytical pre¬ 
dictions. Furthermore, the step decrease of T (t) as a 
function of time in Fig. [TJb) can be explained from an¬ 
other perspective as follows: After a field quench, the 
presence of the external field can not only drive the cor¬ 
related wave packet to move along the bound band but 
also give rise to the Zener tunneling between the bound 
and scattering band of Ho. The Zener tunnelling takes 
place almost exclusively when the momentum of the wave 
packet reaches to the point of center momentum K = 0, 
which is either at the top of the bound band or at the bot¬ 
tom of the scattering band. In this case, the correlated 
pair can be partially pumped into the upper scattering 
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FIG. 7. (Color online) Numerically computed evolution of the functions T (t), r ( t ), and E ( t ) for the initial BP wave packet 
in the form of Eq. with parameter values given in the text and for two typical values of F\ (a) F = —0.097120, and 
(b) F = —0.097815. The other parameter values are U = V = —6.24, k = 1 and N = 111. (a) shows that the evolved state 
remains its correlation and preserves the constancy of the average interval energy. This denotes that the two correlated particles 
undergoes a BO with period ts ~ 65. While in (b), the increase of r(t), E (t) and decrease of T (t) are associated with the 
Zener tunneling from the bound band to the scattering band, which leads to the dissociation of the two-particle correlation. 


band accompanied by the destruction of correlation and 
the increase of energy. At this time, the evolved state 
is a mixed state composed of correlated and scattering 
states. Over a period of time, the scattering and corre¬ 
lated components of the evolved state once again arrive 
at the transition point simultaneously with a character¬ 
istic period given by ts = 27 t/F. Then the two com¬ 
ponents of the evolved state exchange the energy and 
correlation from each other due to the Zener tunneling. 
Through a number of these processes, the T (t) tends to 
be a steady value. The corresponding time tf can be 
termed as quench relaxation time, which deeply depends 
on the distribution of the initial state in the energy space. 
To further understand the field-induced energy variation 
arising from the Zener tunneling from the bound band 
to scattering band, we study the average energy of Hq, 
which can be defined by 

E(t) = (*(t)\H 0 \*(t)). (20) 

As comparison, the average energy E (t) as a function of 
time for two typical quenched systems is plotted Fig. 0 a) 
and (b). In Fig. 0a), the average energy of the evolved 


state oscillates around the strength of on-site interaction 
U = —6.24 preserving the constancy of the average inter¬ 
val energy. The corresponding oscillation period is given 
by tg = 27 t/F. This implies that the two correlated par¬ 
ticles undergo a BO, which is in accordance with our pre¬ 
vious analysis. While for the dynamics of correlated par¬ 
ticles in the quenched system with multi avoided crossing 
regions, the average energy of the evolved state is first 
around the strength of on-site interaction U = —6.24 
and then oscillates around 0, which indicates that the 
evolved state is completely pumped into the scattering 
region with the increase of the average energy of two 
particles. This can be seen from Fig. 0b). On the other 
hand, when the quench field is switched off at time tf. 
the time evolution of the state |T (tf)) is driven by the 
Hamiltonian Hq leading to a constant average energy. At 
this time, the state IT (tf)) are mainly distributed on the 
scattering band of H 0 owing to the increase of the average 
distance r, which indicates that the quench process ac¬ 
complishes the energy transfer from the field to particles. 
This is analogous to the electrothermal effect in the con¬ 
text of classical physics, which describes the interconver¬ 
sion of thermal energy difference and electrical potential. 

















T(F) 



FIG. 8. (Color online) Behavior of the function Tt f ( F ) versus 
the held strength F for the parameter tf = 800. The other 
parameter values are the same with that of Fig. [I] Note 
that the Tt f ( F ) is periodic with period Fb = 0.0015, which 
accords with the period of periodic occurrence of the multi 
avoided crossings in the system. 


Therefore, such a quench process, in which the energy 
of the two particles is elevated by the external field, can 
be termed as a quantum electrothermal effect. Further¬ 
more, one can extend this result to the system consisting 
of multi correlated pairs in which the pair-pair interac¬ 
tion is neglected. The quench field can elevate the total 
energy of the particles associated with the disassociation 
of the pairs in the coordinate space. This process is re¬ 
flected by the characteristic temperature of the system, 
which can be measured in experiments. 

Now we turn to study the sensitivity of the transition 
between the correlated and uncorrelated particles states 
after a field quench. In the previous section, we have 
pointed out that a slight variation of the external field 
can trigger the periodic occurrence of the multi avoided 
crossing regions, which can affect the stability of corre¬ 
lated pair. Thus the dynamics of correlated pair is sensi¬ 
tive to the external field and can characterize the change 
of the energy spectrum. Based on this analysis, we in¬ 
vestigate the time evolution of the initial state |T (0)) 
by suddenly quenching the field strength from an initial 
field-free value to a different value F{t > 0). And we em¬ 
ploy long time-scale transfer rate T tf (F) to measure the 
sensitivity of the correlated-pair dynamics to a quench 
field, which can be described by 


T tf (F) = j2\(rK\e- iH{F)t n*m 


( 21 ) 


where tf is the quench relaxation time. At this time, the 
T tf ( F ) tends to be a steady value. We plot the quantity 
T tf (F) as a function of quench field F in Fig. [H It can 
be observed that the long time-scale transfer rate T tf ( F ) 
from an initial correlated state to final uncorrelated parti¬ 
cle states is sensitive to the quench field strength and ex¬ 
hibits a periodic behavior. And the maximum of T tf (F) 
is corresponding to the initial correlated pair placed on 
the unmixed region of the quenched system. While the 
minimum of T tf (F) denotes the correlated-pair dynam¬ 
ics in the quenched system with multi avoided crossings. 
The numerical results clearly show that the T tf (F) is 
a good indicator of the transition between systems with 
and without multi avoided crossing regions induced by 
the external field. This may also provide a method to 
identify a subtle change of the system parameter in ex¬ 
periments. 


V. SUMMARY 

In this work, we have theoretically investigated the co¬ 
herent motion of two correlated particles driven by an ex¬ 
ternal field in the framework of the one-dimensional Hub¬ 
bard model. We have shown that the particle-particle 
interaction can induce BP states, exhibiting short-range 
quantum correlations. The spectrum comprises two 
Bloch bands, uncorrelated-particle and BP bands. The 
presence of the external field can not only changes the en¬ 
ergy spectrum from the continuous to discrete but also 
bring about the mixture among the correlated state and 
multi uncorrelated states, which is accompanied with the 
exchange of the energy and correlation between the in¬ 
volved two types of states. We have also shown that 
a small change of the external field can induce periodic 
emergence of the multi avoided crossing regions leading 
to the transition from the initial correlated state to fi¬ 
nal uncorrelated states. The corresponding transfer rate 
is sensitive to the quench field strength and exhibits a 
periodic behavior. This process is irreversible and ac¬ 
complishes the energy transfer from the field to parti¬ 
cles, which could be the mechanism of the quantum elec¬ 
trothermal effect. It is envisaged that the sensitivity of 
the correlated-pair dynamics to an external field can be 
applied to detect a subtle variation of system parameter, 
stimulating further theoretical and experimental investi¬ 
gations on the dynamics of strongly interacting particles 
in a variety of driven cold-atom systems. 
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